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r^ ■ Abstract 

O, 
, ^,. An exact calculation of the phase diagram for a loop gas model on the 



brickwork lattice is presented. The model includes a bending energy. In the 
dense limit, where all the lattice sites are occupied, a phase transition occuring 



^Sj . at an asymmetric Lifshitz tricritical point is observed as the temperature 

C""^ . associated with the bending energy is varied. Various critical exponents are 

0\ 



calculated. At lower densities, two lines of transitions (in the Ising universality 
class) are observed, terminated by a tricritical point, where there is a change 
Ch , in the modulation of the correlation function. To each tricritical point an 

o 
o 



associated disorder line is found. 
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I. INTRODUCTION 

Models of closed loops on lattices in two dimensions have attracted considerable attention 
|jl[]. In a theoretical context they arise naturally as high-temperature expansions of spin 
models and they are closely related to integrable systems such as vertex models [0. Loops 
on a lattice may also be regarded as simple models for (short) ring polymers in solution [^Q. 
The segments of the loop are then regarded as monomers, or small clusters of monomers. 
While realistic systems are three-dimensional, the two-dimensional case provides rich critical 
behaviour, and it may be hoped that some features hold in higher dimensions. Solid-on-solid 
models used in the study of roughening transitions in three-dimensional growth may also be 
mapped onto various two-dimensional loop-gas models ||ll. 

In this paper a model of loops is studied, consisting of self-avoiding rings with a bending 
energy. We are mainly interested in the effects of varying the density of monomers on the 
lattice and the temperature. The model is defined as follows. Each bond of the lattice is 
either occupied by a monomer or empty. Each monomer placed on the lattice connects to 
two others such that the only allowed configurations consist of closed self-avoiding loops. 
An energy penalty e is associated with each turn. The density of lattice bonds occupied is 
allowed to change, and the model is studied in the grand canonical ensemble by introducing 
a fugacity K for each monomer. 

On the square lattice, and in the limit that the lattice is maximally occupied (all the 
sites visited) this model corresponds to the so-called F Model [^, which exhibits an infinite 
order phase transition as the bending energy, or equivalently temperature, is varied. At 
low temperatures corners are expelled from the bulk, while at high temperatures there is a 
proliferation of corners. 

A qualitatively similar transition is seen in a model of a single Hamiltonian walk on the 
square lattice with a bending energy [^ . The Hamiltonian walk may be thought of as the 
limit of an interacting self-avoiding walk where the attractive nearest-neighbour monomer- 
monomer interactions are strong enough to exclude any lattice vacancies. In the other limit, 



the non- interacting self-avoiding walk, it is known that the bending energy is irrelevant P,Q ; 
it changes the effective size of a monomer without changing the critical behaviour. At low 
enough densities we expect that for the loop-gas model the bending energy will also be 
irrelevant in this sense. 

One therefore expects that between zero density (bending interaction irrelevant) and 
density one (bending interaction relevant) there should be a "critical" density corresponding 
to a change of behaviour. 

In this paper we study the loop gas model on the brickwork lattice. The brickwork 
lattice corresponds to a square lattice with half the horizontal bonds removed, giving rise to 
the "brick wall" motif, see Figure 1. The brickwork lattice is topologically identical to the 
hexagonal lattice. 

At density one we find a qualitatively similar phase transition to that on the standard 
square lattice in that there is a low temperature phase in which the corners are expelled from 
the bulk and a high temperature phase where there is a finite density of corners in the bulk. 
The details of the transition are however very different: the transition occurs at a tricritical 
Lifschitz point and the high-temperature phase is modulated and critical. Such a phase 
transition is reminiscent of the Pokrovsky-Talapov transition 0]. In this limit our model 
is equivalent to a modified KDP model on the hexagonal lattice PJ^. The existence of a 



modulated phase at p = 1 implies the existence of lines of disorder in the phase diagram [|10 
along which the correlation functions change from modulated at high densities to monotonic 
at low densities. 

As the density is reduced the tricritical point is extended into a line of critical points 
in the Ising universality class which terminates at zero temperature at a critical density of 
about 0.8. For K < 1, or at low density, another line of critical points is observed, again 
in the Ising universality class; our model in this region is essentially the Ising model on the 
hexagonal lattice with one of the three couplings different from the other two. 

The remainder of this paper is organised as follows. In Section ^ the grand canonical 
partition function (and other relevant quantities) is calculated by expressing it first in terms 



of Grassman integrals, which are then exactly computed in the thermodynamic limit. In 
Section |ITI| the i^ ^ oo (p ^ 1) limit is discussed, along with the nature of the low and high 
temperature phases. In Section |IV| the phase boundaries and lines of disorder are found and 
the different aspects of the phase diagram discussed. Section [V] is devoted to final discussions 
and conclusions. 

II. THE MODEL 

We consider a two dimensional gas of loops on a brickwork lattice (BW). The loops are 
self-avoiding, and we assign a fugacity K to each occupied link. This BW lattice can be 
visualized as a honeycomb (HC) lattice (see Fig. 1), and we associate a weight of e~^^ (where 
/? = 1/T is the inverse temperature) to each corner of a loop or equivalently A = e~'^^'^ to 
horizontal links. This model is a straightforward extension of the F-model to the BW lattice. 
The aim of this paper is to calculate the phase diagram and properties of this system as 
a function of the site density p (or equivalently bond fugacity K) and temperature. The 
grand canonical partition function Z of the system is given by: 

Z=J2 K'^Zj, (1) 

Af=0 

= e-^^^(^'^) (2) 

where f2 = L^ is the total number of sites of the lattice of linear dimension L, and Z^q is 
the iV— site partition function. As usual, the canonical partition function at site density 
p = N/Q can be obtained through (|I]) by: 

Z^ = ^^e-^[^i°s^-^^(^'^)l (3) 

In the thermodynamic limit, we have the usual relation 

p = K^\ogZ (5) 



The loop gas can be identified with the graphs of the high-temperature expansion of an 
Ising model on a BW lattice, with a weight K per vertical bond and \^K per horizontal 
bond. This identification holds only provided that K <1. However, the solution follows for 
any value of K. 

Using the results of Houtappel [|ri|,|T^, we have, for any K: 



f{K, T) = -^ r dK r dky \og{a^K* + a^K^ + ag) (6) 

Ibvr^ JO Jo 



where 



al = 1 + 4A^ cos^ ky — 4A^ cos kx cos ky (7) 

a2 = 2 — 4 cos^ ky + 4A^ cos k^ cos ky (8) 

a3 = 1 (9) 

(10) 

This result as well as the correlation functions can be easily obtained through the use of 
Grassman variables. Following Samuel M,IT3, the partition function Z can be represented 



as a Grassman integral 

Z = fVijjVipe-^ (11) 

where 

A \~^ I i m,n i m,n , i m,n i m,n , m,n m,n , m,n m,n /io\ 

^= -2^m i^l +^i^2 +^Z^l + ^4 ^2 (12) 

m,n 

+ A(^1 ' ^2 + ^3 ^4 +'^2 tPa + ^1 ^4 

m,n m,n , m,n m,n , m,n m,n , m,n m,n\ 
V'l ^2 + ¥^3 <^4 + <^2 ¥^3 +^1^4 ) 

I T^/\2 mm I mm , ; m,n+l m,n , mm im— l,n\\ 
+K{X ^3 ' V^i ' + V^4 ' ^2 +^A i'2 )) 

and t/'™'",-?/'™'", ■?/'^'", ■?/'^'" are Fermionic fields attached to each lattice site {m,n). The 
Fermion integral can be performed and the grand potential ^ can be recovered. In addition, 
the generic correlation functions read: 



G(,n,„) = (-ir^l ,kI *„-^^,_^^,_ (13) 

The actual correlation functions contain regular multiplication factors which do not modify 
the long distance behaviour. 

Integration over kx can be performed and gives: 



/( A-, D = -!/'*, log i±^^f^ (14) 

OTT JO I 
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/o 27r Va^ - &2 



G(,„.„) = (-1)" r #^^^ f -'^^^-" ! (15) 



where 



a = (1 + 4A^ cos=^ A;j,)K^ + (2 - 4 cos^ ky)K^ + 1 (16) 

b = -4A=^ cos kyK^{K^ - 1) (17) 

The canonical free energy (at bond density p = N/Q) is given by: 

/,(p,T) =plogK-/ (18) 

where K is determined as a function of p through: 

^ = J_ Tdk 1 dg{K,T,k y) 

K Sir Jo ^g{K,T,ky) dK 

and 



^ = — dky- — - ^^ ' ' ^^ (19) 

K SttJo ^g{K,T,ky) dK ^ ' 



S = ° + 7^ (20) 

In the following, we will work in the grand canonical ensemble, and transpose the results 
to the canonical ensemble when necessary. We first consider the fully-packed lattice (p = 1) 
and then discuss the dilute case. 

III. THE FULLY-PACKED LATTICE 

By analogy with polymer theory, it is interesting to consider the case where all lattice 
sites are visited once and only once by the loops. This is the non-connected version of 
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Hainiltonian path model, with a penalty factor A per corner. From Eq. (|T9D, we see that 
p = 1 for i^ = oo. We are thus led to study Eq. (|18[) in the limit when K ^> go. One 
obtains: 

fc = —-. — 7: dkx I dky\og{l + A\'^ co^ ky — A\^ co^ kx COS ky) (21) 

47rp JO Jo 

As usual we identify the critical points of the system from the zeros of the argument of 

the log in the above equation. It may be seen that no zeros exist for A < l/-\/2 and that for 

A > 1/v^ zeros exist at: 

fc. = (22) 

This implies that the whole region A > l/-\/2 is critical, with a temperature dependent 
critical wavevector. We therefore identify A = l/-\/2 with a tricritical Lifshitz point, and the 
region A > 1/v^ as a Lifshitz line of critical points. Using the definition A = exp — e/T, the 
corresponding temperature for the tricritical point is T^ = 2e/ log 2. 

The integration over k^ in Eq. (pT]) may be carried out explicitly, giving: 



r 1 r^, 1 A+4:X'^COsHy+\l-4:X'^COsHy\. .. 

fc = -^ I dky log( ^— ^) (23) 

or 

1 /-arccos(l/2A2) 

fc = -^l dky\og{AX^cosHy) for T > T, (24) 

= for T<T^ (25) 

(26) 

This form is similar to the models studied by |]^,|^ (see Fig. 2) . 

It is natural to define the average number density of corners n^ as the order parameter 

for this transition. Indeed, we find (see Fig. 3): 

2 1 

nc = - arccos(— — ) for T >Tc (27) 

7T lA^ 

= for T <T^ (28) 



The critical behaviour of the order parameter is given by ric ~ 5T^/^ as T ^ T^ so that 
the critical exponent jS is equal to 1/2. The low-temperature phase is completely frozen, 
consisting of straight vertical lines, with all the corners rejected to the outer boundary. The 
high-temperature phase is modulated in the ^/-direction with a wave vector given by Eq. (^21) 

The same critical behaviour is also seen in the zero-temperature phase diagram of the 
frustrated Ising model on the triangular lattice with appropriately chosen coupling constants 
PJ^ . This model may be mapped onto a tiling consisting of three types of lozenge . One 
lozenge has a lower energy than the other two. At zero temperature, the tiling must be 
perfect. One rapidly realises that the only way of introducing a lozenge of higher energy is 
to introduce an infinite line of them. Identifying the side of a lozenge with the bisector of an 
occupied bond on the dual hexagonal lattice, our loop model may be seen as equivalent to this 
lozenge tiling (see Fig. 4). At non-zero temperature (as defined in our model) a defect line 
may be seen as a restricted SOS interface crossing the lattice. The energy needed to create 
one such line is E = 26L and the entropy is S* = L log 2. When the associated free energy, 
Fi = L{2e — Tlog2), becomes negative, defect lines (and hence corners) proliferate. This 
defines the critical temperature as Tc = 25/ log 2, consistent with the tricritical temperature 
found analytically above. 

In the high-temperature phase, where these lines proliferate, we give a simple physical 
argument for the observed modulation in the correlation functions. The free energy Fi is 
simply the chemical potential for creating one such line. When a finite density of lines 
is present, the reduction of entropy must be taken into account [jT3[, yielding an effective 



repulsion between them. The total free energy for N lines, per occupied bond, is: 

N rp 

F^ = -F,N + cY.Y2^ (29) 

where di is the distance between lines i and i + 1. Minimising F^ with respect to the di, 
subject to the constraint Y^idi = L, gives all the di equal and given by: 



explaining the form of the temperature dependence of the modulation wavevector, Eq. (p2D 
Close to the tricritical Lifshitz point, the free energy scales as: 

f-:;^si^Tr (31) 

from which we obtain the specific-heat critical exponent a = 1/2. 

Along the Lifshitz line, the critical behaviour of the correlation functions can be analysed; 
a generic correlation function is given by: 

Q(^ = __ / / d'^k - (32) 

^ ' 47r2 io io 1 + 4A4 cos2 ky - 4A2 cos k^ cos ky ^ ' 

up to a regular multiplication factor. 

Away from the tricritical point, i.e. A > l/"\/2, we may develop around the critical wave 
vector (0, k) defined in Eq. (0) Setting k^ = qx and ky = k + qy, this becomes: 

The exponential prefactor gives the expected spatial modulation in the y-direction, and the 
correlation function has a logarithmic behaviour at large distances. 

Around the tricritical point, the critical wave vector vanishes as well as the coefficient of 
the q^ term, and the expansion must be carried to the next order: 

where /i = 2A^ — 1. Therefore, the correlation functions have anisotropic scaling, with critical 
exponents u^ = 1 and Uy = 1/2. 

IV. THE DILUTE LATTICE 

We now move to the dilute case p < 1. 



A. Critical lines 

As the density p or the fugacity K is lowered, the tricritical Lifshitz point extends into 
a critical line. This line can be obtained from the zeros of the logarithm of Eq. (^ 

A;^ = 

K = , -^ (35) 

VI - 2A2 ^ ' 

This critical line exists for K > 1. In the p — T plane, its equation close to the fully-packed 
case p = 1 is given by: 

P ^ 1 + 2e^ (36) 

c 

where T^ = 2e/log2 is the Lifshitz tricritical temperature. 
For fC < 1, there exists another critical line given by: 

fc^ = 

K = ^ ^ (37) 

VI + 2A2 ^ ^ 

In this region (K < 1), the fugacity can be identified with the tanh/3Ji of a regular 
anisotropic Ising model on a honeycomb lattice. Similarly, the second coupling is given 
by X^K = tanh/3J2. 

It can be easily seen that both lines correspond to the 2d Ising universality class : z/ = 1. 
Note that since the problem is formulated as a loop gas, the correlation functions don't 
correspond to the spin correlation functions of the Ising model; here we have r] = 0. 

The phase diagram in the K — T plane is shown in Fig 5. Using equation |1^, we find the 
phase diagram in the p — T plane (see Fig. 6). One can identify a high-density transition 
line and a low-density transition line. The high-density line ends at pc — 0.8185 and the 
low-density one at p^ ~ 0.19. 
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The two critical lines merge at K = 1, where three phases become critical simultaneously, 
defining a tricritical point. This is manifested in the p — T plane by a jump in the critical 
density p at K = 1. As usual for zero temperature tricritical points, observables develop 
essential singularities. 

B. Disorder line 

/^From Eq. ([15D , it is easily seen that the correlation functions change from oscillating 
(in the x-direction) for K > 1 to monotonic for K < 1. These two regimes must therefore 
be separated by a line where the short distance correlation changes from oscillating to non- 
oscillating. This line K = 1, which passes through the tricritical point {K = 1,T = 0), is 
called a disorder line [0. According to the definition of Garel and Maillard |T^, this is a 
line of disorder points of the first kind (with zero correlation length). 

We have seen that at p = 1, there is a Lifshitz critical point separating a frozen low- 
temperature phase from a modulated (in the y-direction) high-temperature phase. At lower 
densities, the correlation functions are not modulated. This happens separately for each 



value of kx- Following Garel and Maillard ||T6[, we define the disorder line as the line for 



which the first mode [k^ = 0) changes behaviour: 

K = ^ ^ (38) 

V2A2 - 1 ^ ^ 

This line is defined in the high-temperature region A > l/-\/2 only, and corresponds to a line 
of disorder points of the second kind. 

V. CONCLUSION 

In this article, a loop gas model on a brickwork lattice was considered. An energetic 
penalty was included for each corner. At p = 1 we observed a phase transition from a 
low-temperature frozen (corner free) phase to a high-temperature phase modulated in the y- 
direction. The phase transition occurs at a tricritical Lifshitz point, where z/^; = 1, z/^, = 1/2. 
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The whole high-temperature phase is criticaL These results are reminiscent of a phase 
transition of the Pokrovsky-Talapov type. This behaviour is completely different from the 
critical behaviour of the analogously defined model on the square lattice at p = 1 (the F- 
Model). This is due to the combination of two effects; the brickwork lattice automatically 
imposes self-avoidance without the inclusion of additional fugacities, and the brickwork 
lattice is intrinsically asymmetric. 

The phase diagram is given in the K — T (and equivalently the p — T) plane. Two lines of 
critical points were observed corresponding to high and low-density phase transitions. The 
high-density phase transition is to a phase modulated in the x-direction, and the low-density 
phase corresponds to the usual Ising transition. Both transitions are in the Ising universality 
class and meet at T = at another tricritical point. 

While the model studied is simple, the resulting phase diagram is suprisingly complex. 
In the formalism chosen, it is not clear how to characterise the different finite-density phase 
transitions in terms of the loop-model observables. 
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FIGURES 



FIG. 1. The Brickwork and the corresponding Honeycomb lattice. 
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FIG. 2. The free energy as a function of temperature for the fully packed model. 
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FIG. 3. The average number of corners as aJunction of temperature for the fuhy packed 



model. 
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FIG. 4. The mapping between the loop gas and the lozenge model: the shaded region shows 
an infinite line of defects. 
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FIG. 5. Phase diagram of the loop gas in the K — T plane; the solid lines denote the phase 
boundaries; the thick line represents the Lifshitz line; the dashed lines correspond to the disorder 
lines. The three phases are defined in the text. 
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FIG. 6. Phase diagram of the loop gas in the p — T plane; the solid lines denote the phase 
boundaries; the thick line represents the Lifshitz line; the dashed lines correspond to the disorder 
lines. The three phases are defined in the text. 



19 



